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The weak coupling region of CP^ -1 lattice field theory with the #-term is investigated. 
Both the usual real theta method and the imaginary theta method are studied. The latter 
was first proposed by Bhanot and David. Azcoiti et al. proposed an inversion approach 
based on the imaginary theta method. The role of the inversion approach is investigated in 
this paper. A wide range of values of h — — Im6 is studied, where 9 denotes the magnitude 
of the topological term. Step-like behavior in the x-h relation (where x = Q/V, Q is the 
' topological charge, and V is the two dimensional volume) is found in the weak coupling 

qq | region. The physical meaning of the position of the step-like behavior is discussed. The 

, inversion approach is applied to weak coupling regions. 

o 

q ' §1. Introduction 

The two-dimensional lattice CP^ 1 model with a #-term is investigated. The 

problem of obtaining the partition function Z{6) numerically stems from the difficulty 

in treating the complex valued Boltzmann weight. This difficulty is avoided by 

Q-i, expressing Z{6) as a Fourier seriesP 
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z{e) = Y J P{Q> 



where P(Q) is the topological charge distribution, i.e., the probability of finding a 
topological charge Q in the system at 6 = 0. 

In the strong coupling region, P(Q) can be approximately expressed as a Gaus- 
sian function, P(Q) oc exp(— yQ 2 ), and a first-order phase transition at 6 = ir 
is obtainedP - ^ In the weak coupling region, P(Q) exhibits behavior that differs 
greately from the Gaussian form. Instead of a quadratic Q dependence, an almost 
linear form is founcP'^ in the exponent of P(Q): 

P(Q) occ IQI =e IQ|lnc . (1-1) 

In the weak coupling region, c was found to be a quite small constant. This linear 
exponent is a simplified typical form. 
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***' E-mail: yoney ama@cc . saga- u . ac . j p 
f) See Eq. (4.9) of Ref. 0. 
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Bhanot and David first proposed the imaginary theta method, ' in which the 
parameter is taken to be purely imaginary, 

6 = -ih, 

with h being a real parameter. Azcoiti et al. proposed an inversion approach 71 based 
on the imaginary theta method. 

In the imaginary 6 case, numerical simulations can be performed for Z(h), since 
the Boltzmann factor becomes real in this case, and we have 

where S is an action and z (z*) denotes appropriate fields (their complex conjugates). 
But it seems that the meaning of the inversion approach 7 ' is not well understood. 
For this reason, the role of the inversion approactP based on the imaginary theta 
methocP is investigated in this paper. 

We performed a numerical analysis with 9 = —ih for both strong and weak 
coupling regions. After presenting the results of this analysis, we discuss the meaning 
of the approach used by Azcoiti et al. Our understanding of the meaning of their 
approach is summarized as follows. 

1. In some cases, real theta results can be obtained from imaginary theta results by 
analytic continuation. However, this is not true in other cases. What, then, does 
the imaginary theta method mean? It does not mean analytic continuation at 
nonzero theta. The inversion approach is nothing but one of the fitting methods 
of the topological charge distribution P(Q) at 8 = . This is shown in §3. 

2. The imaginary theta method is suited to determining the ^-dependence for a 
wide range of values of h. In particular, the x-h relation for a wide range of 
values of h, and thus a wide range of values of x, is obtained, where x = Q/V. 
From this x(h)-h relation, the inversion approach leads to the h{x)-x relation. 

3. In the strong coupling region, the Gaussian form of P{Q) is reconfirmed using 
the inversion approach. 

4. In the weak coupling region, we have found "step-like behavior" in the x-h re- 
lation. The position of a step gives the value of the parameter c(f3), where (5 
represents the inverse coupling constant of the CP^ -1 model. The parameter 
c((5) is that found in our previous analysis of the topological charge distribu- 
tiorP at = 0, P{Q) oc c(/3) Wl. 

5. In the weak coupling region, fluctuation of the field variables z and z* is greately 
suppressed, and thus the probability of topological charge excitation is also 
greately suppressed in regions of small h. The parameter c{(5) is a measure of 
the degree of suppression of the topological charge excitation. 

6. The h{x)-x relation is obtained from the x(h)-h relation in the inversion ap- 
proach. The functional form of h(x) is obtained by fitting to the data with 
the appropriate function. Actually, at the end of f(x) in the weak coupling 
region is obtained by integrating h{x) = df{x)/dx in the inversion approach 
[see Eq. (|3-22[) ]. and the leading term in that result reproduces the result of the 
linear model [exponent of Eq. ]. 
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This paper is organized as follows. The inversion approach based on the imag- 
inary theta method is explained in §2. The results of the numerical calculation are 
presented in §3. Conclusions and discussion are given in §4. 

§2. Inversion approach in the imaginary theta method 

2.1. Formulation 

Much progress has been made in non-Abelian lattice gauge theory; e.g., asymp- 
totic freedom was confirmed by the observation of string tension.^ It appears possible 
that quark confinement can be explained as an "area law" in the lattice formulation. 
Instanton excitation allows the existence of the topological term, namely a theta 
term in the action. However, lattice field theory with the theta term is not well 
understood, because the Euclidean formulation introduces a complex Boltzmann 
factor, and it does not allow direct Monte Carlo simulation. Bhanot and David 
first introduced a purely imaginary theta parameter and studied the 0(3) non-linear 
sigma modelP If we take 9 to be purely imaginary, the Boltzmann weight becomes 
a real positive quantity, and this allows dirfect numerical simulation. Azcoiti et al. 
introduced the inversion approach based on the imaginary theta method in Ref. |7J). 
However it seems that the meaning of the inversion approach is not well understood. 
In order to understand the meaning of the inversion approach, we employ the imag- 
inary theta method in both the strong and the weak coupling regions and study the 
inversion approach. We then compare the imaginary and real theta methods and 
study the role of the inversion approach. 

We begin with the real 9 case. The CP^ -1 model with the #-term on a two- 
dimensional Euclidean lattice is considered!*^ The action with the #-term is defined 
by 

S e (z,z*) = S(z,z*)-i9Q(z,z*), 

where 

S(z, z*) = { 1 " E K(n>a(n + fi)\ 2 \ (2-2) 

is the action and z a (n)(a = 1, • • • N) denotes a CP^ -1 field on each site n. The site 
n + pi, is the site nearest to n in the direction fi. The topological charge Q(z,z*) is 
defined as 

Q( z , z *) = -Lj2 A n, (2-3) 

A D {z, z*) = - Y,{ A M) + A ^ n + A) - M n + Mn)}e^, (2-4) 

where the quantities An{n) are defined as 

exp(i74^(n)) = z\n)z{n + jj,)/\z\n)z(n + £l)\ (2'5) 



2-1 
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by the CP^ -1 field z. For the CP^ -1 model, the complex field z a (n) satisfies the 
equation 

N 

z ] (n)z(n) = z* a (n)z a (n) = 1. (2-6) 

a=l 

The partition function for the two-dimensional CP^ -1 field theory is given by 

J VzVz* exp(-S(z,z*) +i9Q(z,z*)) 

J VzVz* exp(—S(z, z*)) 

J (VzVz*)^ exp(-S(z, z*) + i0Q(z, z*)) ( 2 . 7 ) 
= ^ 1 

Q=integer / VzVz* exp( — S(z, Z* )) 

= Yl P (Q)e l9Q , 

Q=integer 

where (VzVz*)^ denotes the constrained measure in which Q(z,z*) = Q. Here, 
P(Q) is the topological charge distribution estimated with the action for 9 = as 



Z V {0) 



(VzVz*) iQ) exp(-S) 

P(Q) = —r • (2-8) 

/ VzVz*exp(-S) 
It satisfies the relation 

£i>(Q) = l. (2-9) 

Q 

Once the topological charge distribution is known, the partition function at any 9 
can be obtained from the Fourier series for the case of real 9: 

Z V (9) = Y J P{QY dQ =Y,exp(-Vf v (x Q )y 6V ^, (2-10) 

Q X Q 



where 



XQ = Q/V, f v (x Q ) = ~lnP(Q). (2-11) 



Now we introduce an imaginary 9™ Setting 9 = —ih (withh real), we have 

/ VzVz* exp(-S(z, z*) + hQz, z*)) 



Z v (h) 



^ fVzVz* exp(-S) (2 12) 

^2exp(-VMx Q ))e hV ^ 

X Q 
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Thus h plays the role of the external source for the topological charge xq. Once a 
constant background field h is given, the expectation value x(h) of the topological 
charge per unit volume is given by 



x{h) 



v:D Z *^l eM -s h ) x; XQ P h (Q) 



V 



VzVz* exp(-Sh) 



Yjxq P h(Q) 



(2-13) 



where 



and 



P h {Q) = J (VzVz*)^eMS h )/ J VzVz*eMS), (2-14) 



S h (z,z*) = S(z,z*)-hQ(z,z*). 



(2-15) 



From this form, we can conclude that the imaginary theta method is a kind of 
"trial function" (subtraction) method, in which the action is replaced by that with 
the subtraction term Strial 

S'eff = Sfi = S — Serial- 

The "trial (subtraction) function" is taken as a special form in the imaginary theta 
method, 



-hQ(z,z*). 



Now we explain the inversion approach. When the volume V is large, xq is 
almost continuous, and xq in the sum in Eq. (|2- 12|) can be approximated by xq, 
the value of xq at which exp{— V(/v(xq) — hxQ)} is maximal. This saddle point 
method gives 



with 



where 



Z v {h) oc exp(-V(fv(xQ) - hx Q )), 
dg(xQ) 



dx 



Q 



at xq = xq, 



Equation (|2-17|) gives 



9(xq) = fv(x Q ) - hx Q . 
dfv(xq) 



dx c 



h. 



(2-16) 



(2-17) 



(2-18) 



(2-19) 



The quantity xq is the expectation value of the topological charge (per unit volume). 
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1) The expectation value xq for a given background h can be obtained by numer- 
ical simulation. Specifically, x is regarded as a function of h, and the x(h)-h 
relation can be obtained. 

2) On one hand, due to Eq. Q2-19|l . h is the first derivative of fv(x) at xq. Here- 
after, the suffix V is omitted. 

3) We plot an illustrative example of function of h in Fig. 1. Exchanging 
h and x, we obtain Fig. 2. This is the "inversion" of the x(h)-h relation to the 
h(x)-x relation. In other words, h is now a function of x. 




Fig. 2. Exchanging x and h in Fig.l, we obtain 
. 1. An illustration of the expectation value , rir \ r c 

1 h = J (x) as a lunction ol x. 

X cLS £L function of the background field h. 



From the relation 

df(x) 



h = , (2-20) 

dx x=x 

we have df(x)/dx as a function of x. That is, fitting h by an appropriate function 
of x, we obtain the functional form of df(x)/dx. By integrating this over x, 

df(x') 



dx' 



-dx' = f{x), (2-21) 



we find f(x). In this way, we obtain the topological charge distribution at 9 = 0. 

We should note that this inversion approach is applicable only when x(h) is a 
monotonic function of h. Otherwise, the inverted function h(x) is multivalued. We 
do not treat this case in the present paper. 

At 6 = 0, the function f(x) is usually directly evaluated through numerical 
simulation of the topological charge distribution. By contrast, f(x) is determined in 
the method proposed by Azcoiti et al. as an implicit function through the process 
represented by ()2T6|) through (|2-21j) . Thus, the x-h relation is given by ()2T9|) . 
This implies that the functional form of h(x) is obtained numerically. The relation 
h{x) = df(x)/dx leads to f{x) after integration over x. 



9 -term, CP 1 Model and the Inversion Approach in the Imaginary 9 Method 7 

2.2. Qualitative difference between strong and weak coupling behavior 

In a previous paperP P{Q) of the two-dimensional CP 2 model is numerically 
obtained. Now we study the results from the point of view of the imaginary theta 
approach. In the strong coupling region {(5 = small < 1), P{Q) is approximately 
given by the Gaussian form 




0.1 0.2 0.3 0.4 0.5 0.6 

X 

Fig. 3. Strong coupling case lasa function of Fig. 4. Exchange of x and h in Fig. 3. 

h. 



P(Q)=e~vQ 2 = e - Vax Q, (2-22) 

and 

f(x) = ax 2 . (2-23) 

In this case, we have 

h = f( x ) = lax. (2-24) 

Hence, a linear relation between (h and x) is expected. 

In the weak coupling region, the x-h relation is expected to exhibit "step-like 
behavior" (Fig. 5). In order to make it easy to understand why Fig. 5 is expected, 
we consider typical simplified behavior of the Q dependence in the weak coupling 
region,^ with the exponent of P(Q) assumed to be proportional to \Q\: 



P{Q) ~ C ' Q ' = e' Q ' lnc = e — 1*31 ln (i/c) _ e -v\x Q \\n(i/c) ^ (2-25) 

The parameter c is known to be quite small from phenomenological considerations. 
Thus, in this case, we have f(x) = \x\ m(l/c),^ which we call the "linear x model" 
hereafter. Then the relation h = f'(x) is given by 

h=f(x) = \ ; n j;/ c ! =\ x *\ (2.26) 

J w 1 -ln(l/c) = -h , x<0. y ' 



8 M. Imachi, H. Kambayashi, Y. Shinno and H. Yoneyama 



The behavior in the x > region is shown in Fig. 6. Then, interchanging x and 
h, we obtain the h-x relation, shown in Fig. 5. From Eq. ()2-26j) . Fig. 6 is obtained 
first. But from the actual numerical simulation, x as a function of the background 
h is computed first. Thus, Fig. 5 is put before Fig. 6. 

We set fes( x ) = f{x) — hx. By employing the simplified functional form "linear 
x model", f{x) = \x\ho, we investigate the three cases (i) — (iii) below. 

(i) < h < ho case: 
We have 

= ( 7u k " M % ' X< n 0, ( 2 - 27 ) 
\ {ho - h)x, x > 0. v ' 

Then exp(— Vf e s{x)) is peaked at x = 0, and we expect x ~ 0. 

(ii) h = ho > case: 

In this case, we have 

/cff(x) = \ 0, x>0. (2 ' 28) 

Then exp(— Vf e g{x)) is constant in the x > region, and we expect x will 
take a positive value between and the maximum possible value; that is, x is 
undetermined. 

(iii) h > ho > case: 
Here we have 

J K ' \ -{h-h )x, x > 0, v ; 

where h — ho is positive. Then exp(— Vf e g(x)) favors as large a value of x as 
possible and we expect that x is given approximately by the maximum possible 
value. Since Q is bounded from above in the finite volume case, we have x < 
Qma,x/V. With periodic boundary conditions, \Q\ < ^ is the limitP Then \xq\ 
is bounded from above by xo (< 1/2). 

Summarizing (i) — (iii), we have found the qualitative behavior of x in the weak 
coupling region as shown in Fig. 5, namely, "step-like behavior" of x as a function of 
the background h. 

The expectation value x is expected to be quite small in < h < ho region, 
while x becomes approximately xq in the h > ho region. The position of the step- 
like increase of x is expected to be located at h = ho = ln(l/e). 

We now give a short comment. In the analysis of finite density QCD, the chem- 
ical potential fi and the nucleon density n (= N/V) (where N and V denote the 
nucleon number and the volume of the system) enter in place of the parameter h 
and the topological charge density x. The low temperature, T ~ 0, regime of QCD 
corresponds to the weak coupling, large /3, regime in the CP^ -1 model, where the 
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2 4 6 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



X 



Fig. 5. Weak coupling case. The x-h relation Fig. 6. A simplified model in which fix) = 
is inferred from Fig. 6. A step-like increase |x|ln(l/c) gives h = ho = ln(l/c) ~ 10.0 

at h = ho is expected. The value ho is for x > 0. 

tentatively taken to be 10.0 here. 



topological charge excitation is suppressed due to the fact that only a small fluctua- 
tion of the gauge field is allowed. The step-like behavior schematically shown in Fig. 
5 is expected to occur in finite density QCD analysis in T ~ oP"'*) 

§3. Numerical calculation of the CP 2 model 

3.1. Numerical results 

(1) As described in §2, the expectation value x(h) was obtained by numerical sim- 
ulation. The expectation values x(h) were calculated at various values of the 
background source points h. 

(2) The saddle point method gives the relation 

h = f'(x). (3-1) 

The above (1) and (2) give the relation between x(h) and h at many points. From this 
x-h relation, the form of f'(x) as a function of x is obtained by fitting the calculated 
points. Once the functional form /'(£) is obtained from this fitting process, we obtain 
f(x) itself by integrating /'(£), 



f{x) = / (3-2) 

J 

*' The step-like increase of n at fi = fio is schematically given. [See Figure 8.10(a) in the textbook 
of Kogut and Stephanov.] See also the paper by S. Kratochvila and Ph. de Forcrand^^^ The 
results (i) — (iii) of this section are quite similar to the results presented in Fig. 2 (the T < Tq 
case) of Ref . 11 U . 
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In Ref. HJ), we discussed the "direct method" and the "indirect method". The 
latter is equivalent to the fitting method. The imaginary theta method is simply a 
candidate for the "indirect method". 

Now we present the results of the numerical simulation of the CP 2 model using 
the imaginary theta method. For each (3 and h, the number of measurements was 
set to 10 5 . 



. 45 



X 




18 20 



h 

Fig. 7. Expectation value x for various values of the background h. The values of the inverse 
coupling j3 are 0.0, 1.0, 2.0, 3.0, 4.0 and 5.0. The lattice size is L = 50 in all cases. 

Figure displays the expectation value x with the volume V = L 2 (L = 50) 
for various values of h = — Im(9 at j3 = 0.0,1.0,2.0,3.0,4.0 and 5.0. In the strong 
coupling cases {(3 = 0.0, 1.0), x increases linearly as a function of h from the origin: 



1 

2a 



h. (fc = 0.0~5.0) 



(3-3) 



For much larger values of h, x exhibits convergence due to the restriction that 
the topological charge cannot exceed V/2 in a finite volumeP Thus x is bounded 
from above: 



(3-4) 



Step- like behavior is observed in the weak coupling regions. In weak coupling 



cases (j3 > 5.0), x is strongly suppressed in comparison with the strong coupling 
cases in the region of small h (h < 5.0), 



(3-5) 
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In regions of larger h, x begins to increase rapidly and reaches x ~ xq/2 at 
h ~ ho(P). In h > ho(/3) region, x begins to converge to a constant due to the 
restriction x < xq. 

The position /io of the step-like increase depends on the coupling constant; that 
is, /io is a function of (3. The observed step-like behavior is not as sharp as the 
simple sharp step-like increase mentioned in §2 (Fig. 5), but it is clearly observed. 
For (3 = 5.0, for example, we have 




Fig. 9. Expectation value x for various values 
Fie. 8. Expectation value x in the strong cou- r , . , ,. , n - „n 

° ot h m the weak coupling case (p = 5.0). 

pling case (B — 0.0). Lattice sizes are r „. . T 1A „„ „„ ,„ , rr , 

r ° y ' Lattice sizes are L =10, 20, 30, 40 and 50. 

L =10, 20, 30, 40 and 50. 



Tin = hi(l/c) ~ 10.0. (3-6) 

Note that all x in Fig. 7 are monotonic functions of h. 

Figure 8 displays x in the strong coupling case (/3 = 0) for various sizes, L = 
10, 20, 30, 40 and 50. The values of x for different sizes coincide and exhibit linear 
dependence on h. Because h = f'(x), a linear dependence of x on h for h < 5.0, i.e., 

x = —h, (3-7) 



gives 



Integrating this, we obtain 



2a 



f'(x) = lax. (3-8) 



f(x) = ax 2 , (3-9) 
namely a Gaussian distribution for P{Q): 



P(Q) a exp(-V7(x)) = exp(--Q 2 ). 



(3-10) 
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Figure 9 displays x in the weak coupling case (/? = 5.0) for various sizes, L = 
10, 20, 30, 40 and 50. Smooth step-like behavior is found for all these sizes. More 
detailed behavior for small h is shown in Fig. 10. For small h, some dependence 
on the value of L is observed. The global behavior for the values h = 0.0 - 20.0, 
however, is almost the same for all values of L considered here(Fig. 9). 




Fig. 10. Expectation value of x in the weak 
coupling case (/3 = 5.0) at smaller h (0.0 ~ 
5.0). The lattice size dependence is clear 
for smaller L (10 - 30). The lattice size 
dependence is weaker for L — 40 and 50. 



Fig. 11. Expectation value of h = fix) as a 
function of x for weak coupling (/3 = 6.0) 
at L — 50. Both experimental data and a 
fit [Eq. I|8-21fl ] are shown. 



In Fig. 11, h vs x in the case of weak coupling, (3 = 6.0, for L = 50 is shown. 
Plateau-like behavior of h at h = ho ~ 11.0 is clearly seen. It is smoother than that 
in the simple case considered in the previous section, but Fig. 11 is the reminiscent 
of the plateau-like behavior (Fig. 6). To see the detailed behavior for values of h 
near the origin, a log-log plot is shown in Fig. 12. 



3.2. Real 9, imaginary 9 and analytic continuation 

Employing a simple "linear x model", we now investigate the role of the imagi- 
nary theta method. The linear x model is defined by 



f(x) = \x\ho, ho = ln(l/c) > 0. 



Hence, the exponent f{x) of the topological charge distribution is a linear function 
of the topological charge (divided by the volume). In the case of real 9, we have 

Z v (9) = Y^e- y f^e ivex - 

— e -V\x n \ho e iV0x n _ (3-11) 
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In the case of imaginary 9,8 = —ih is inserted, and we have 
Z v (h) = J2e- vlXnlho e Vxnh 

= (negative Q) + (Q = 0) + (positive Q) 

N N 
_ e -hon-hn _|_ ]_ _|_ e -hon+hn 

n=l n=l 
N N 

= j> n +i+xy, (3-i2) 

n=l n=l 

where x n = n/V = Q/V , s = e~ ho+h and t = e~~ h °~ h , and the finite volume imposes 
an upper bound of the topological charge Q at N. 
Then the partition function is given by 

Z V (h) = ±— t + 1 + 1 S _ s ■ (3-13) 

In the case of real 9, s(9) = e - ho e ie and t{6) = e~ h °e- ie lead to \s(9)\ < 1 and 
\t(9)\ < 1. Thus, s N+1 and t N+1 approach zero, and Eq. ()3-13|) becomes 



l-t(5) 

(3-14) 



1-c 2 



1 -2c cos (9 + c 2 

where c = e _fi ' 1. 

In order to address the question of whether it is possible to extend real 9 to 
imaginary 9, let us study the following two cases. 
(1) Imaginary 9: 9 = —ih, with \h\ < ho 



In this case, s = e ho+h < 1 and t = e h ° h < 1, and both s N+1 and t 1 ' 1 ' can 



— o \ J- , aim uwui o anu 

be safely ignored. Then we obtain 



Zv(h) -> + 1 + S 



1 -t l-s 

1-c 2 . > 

= 9- (3T5) 

1 — 2c cosh n + cr 

Equation ()3T5|) is then simply analytically continued to Ea. ()3T4|) by taking 
h -» *0. 

(2) Imaginary 9: 9 = —ih, with \h\ > ho 

We consider h > ho > 0. It should be noted that 



t = e~ h °- h < 1, s = e - ho+h > 1. (3-16) 

In this case, s is greater than unity. Thus, the leading contribution to Zy{h) is 
given by 
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_ S N+1 



Z v (h) ~ . (3-17) 

1 — s 

From Eq. (|3-17|) . the expectation value of x is 

1 dZ/dh N „ , , 

*=VZ = finite = x < 1/2. (3-18) 

An important point is that Zy(h) is not simply obtained in the h > ho region 
by analytic continuation from Eq. (|3-14f) by taking 9 — > —ih. Rather, Zy{h) is 
given by Eq. (|3T7|) . which is completely different from the analytic continuation 
form Eq. ()3T5j) in the h > ha region. 



This simple "linear x model" provides the important lesson that there is a case in 
which imaginary 9 does not yield the real 6 result by analytic continuation. Rather, 
the imaginary 9 method is used as a fitting procedure for the topological charge 
distribution P(Q) at 9 = oP 

In the weak coupling region, x exhibits step-like behavior as a function of h. 
For h < ho = ln(l/c), x is close to zero. Then, it increases suddenly at h ~ ho and 
then satisfies x ~ xq for h > ho. Actually, the h-x relation obtained numerically 
does not exhibit an abrupt step- like increase at h ~ ho, but a somewhat gentle one. 
Simple functional form representing this "gentle" step-like behavior is expressed as 
a function of h as 



Xoe c e (h-h ) 
x == 

I _|_ e c e (h-h ) ' 

where x is written The parameter values are 

h ~ 11.0, x ~ 0.4031, c e ~ 0.95 

for p = 6.0 and L = 50. 

This relation is easily inverted. We have 



(3T9) 



e c e (h~h ) = ^ 

and 



e c e {h~h ) = (3 . 20 ) 

Xq — X 



h = ho + —{lnx-ln(xo-x)}. (3-21) 
Since h is given by h = f'(x), we integrate Eq. (j3-21|) and obtain 

f(x) = [ f{x')dx' = [ h(x')dx' 
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Fig. 13. Integrating Eq. 13-21H . we obtain the 
Fig. 12. In order to see the detailed behavior analytical form f(x) [Eq. (£23] for this 

of Fig. 11 for smaller h, a log- log plot is cage 
shown for the same set of numerical data 
as in Fig. 11. 



= Hqx H {xlnx + [xq — x) ln(xo — x)} + d, (3-22) 

where d is an integration constant. If we assume /(0) = 0, d becomes 

d = — xolnxo. (3-23) 

c e 

The value of h for (3 = 6.0 and L = 50 is plotted as a function of x = x in Fig. 11 
and 12. In these figures, the fitting function given in Eq. l)3-2ip is also shown. The 
integrated f(x) is plotted as a function of x in Fig. 13. An almost linear dependence 
of f(x) on x is observed in the weak coupling region. 

§4. Conclusions and discussion 

The inversion approach^ based on the imaginary theta method^ was investi- 
gated in both strong and weak coupling regions. In the strong coupling region, the 
expectation value x exhibits linear dependence on h for h < 5.0. In the weak coupling 
region, x is much smaller than that in the strong coupling region: 

ic(weak) <C x(strong) for h < 5.0. (4T) 

The expectation value x is expected to be quite small in the region < h < ho, while 
x is approximately equal to xq in the region h > h$. The position of the step-like 
increase of x is expected to be located at h = ho = ln(l/c). 

The expectation value x thus displays a step-like increase at h ~ ho = ln(l/c) ~ 
10.0 ~ 11.0. The parameter c represents the probability of a single topological charge 
excitation in the weak coupling region. The value of c is quite small. 

We have numerically calculated x for each h, thus obtaining function of 

h. Then we have obtained h as a function of x by inverting x and h. In the large 
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V limit, h is identified with f'(x), and x is simply written as x. If we fit h with an 
appropriate functional form hst(x), then f(x) is obtained by integrating hfa(x) as 

/(*) = J X f'(x')dx' = J* Kt(x')dx' , (4-2) 

where f(x) denotes f(x) = — V~ l lnP(Q). In this way f(x), namely, P(Q) at 6 = 0, 
is obtained. 

This process shows that the inversion approach in the imaginary theta method 
is not the analytic continuation from h to non-zero theta, but P(Q) at 9 = is 
obtained as one of the products of this approach. A simplified fitting function h^ t (x) 
is presented in §3, and the result, f(x), is given for that fitting function. The obtained 
function f(x) [see Eq. (|4-2jl ] is shown in Fig. 13. 

The purpose of the present paper is to clarify the meaning of the inversion 
approach proposed by Azcoiti et al. For this purpose, we have chosen a simple model, 
the CP 2 model. In our previous analysis,^ it was shown that this model exhibits 
qualitatively different behavior of P(Q) in the strong and weak coupling regions. This 
difference emerges as that in the x-h relation of the imaginary theta method. Since we 
have employed the standard action, and this model is contaminated by dislocations, 
precise information about continuum physics cannot be obtained. Although a further 
investigation of the continuum limit is left for a future study, what we have clarified 
here, i.e., that the x-h relation in the weak coupling region exhibits step-like behavior, 
should not be altered if a more realistic model were employed. For example, in our 
previous analysis^ of the CP 3 model with a fixed-point action, it is shown that the 
topological quantities, such as P(Q) and the expectation value {Q)q, exhibit nice 
scaling behavior, and that P(Q) behaves differently in the strong and weak coupling 
regions. Specifically, we studied the behavior of the effective exponent of P(Q), "f e fi, 
and we found that 7 e g in the weak coupling region is smaller than that in the strong 
coupling region (Gaussian). This suggests the possibility of step-like behavior of x-h 
in the continuum limit. Further study of this point is needed. 

"The step-like increase" of x at h ~ ho in the weak coupling region (1/(3 ~ 0) 
is shematically shown in Fig. 5. The results of the actual numerical simulation are 
shown in Fig. 9. As stated at the end of §2, the h-x relation is quite similar to the 
fi-n (chemical potential vs nucleon density) relation in QCD. For T ~ (where T 
denotes temperature), a step-like increase at \i ~ (jlq is expected, as shown in Figure 
8-10(a) of the textbook of Kogut and StephanovP^ Further investigation of the 
correspondence between h-x in the CP^ -1 and fi-n relation in QCD is an interesting 
problem. 
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